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The structure of the large scale distribution of the galaxies have been widely studied since the 
publication of the first catalogs. Since large redshift samples are available, their analyses seem 
to show fractal correlations up to the observational limits. The value of the fractal dimension(s) 
calculated by different authors have become the object of a large debate, as have been the value 
of the expected transition from fractality to a possible large scale homogeneity. Moreover, some 
authors have proposed that different scaling regimes might be discerned at different lenght scales. 
To go further on into this issue, we have applied the correlation integral method to the wider sample 
currently available. We therefore obtain a fractal dimension of the galaxy distribution which seems 
to vary by steps whose width might be related to the organization hierarchy observed for the galaxies. 
This result could explain some of the previous results obtained by other authors from the analyses 
of less complete catalogs and maybe reconcile their apparent discrepancy. However, the method 
applied here needs to be further checked, since it produces odd fluctuations at each transition scale, 
which need to be thoroughly explained. 
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INTRODUCTION 



Standard cosmology is based on the assumption 
that the Universe is spatially homogeneous, at least 
on scales sufficiently large to justify its approximation 
by a FLRW model. 

However, the consensus on a homogeneous feature 
of structures, even on very large scales, has never 
been complete. At small separations, data worked 
out using the correlation function method show a 
correlation lenght r = 5ft- -1 Mpc and the galaxy dis- 
tribution exhibits a fractal structure with dimension 
D 2 ~ 1.2 [l|. The correlation integral method |2J 
gives a fractal dimension slightly larger, D 2 ~ 1.3 
to 1.5. At larger scales, the value D 2 ~ 2 has been 
proposed up to scales of at least 150ft-~ 1 Mpc 0. 
For review articles see Martinez, 1999 [Bj or Wu, 
Lahav and Rees, 1999 0. It has also been suggested 
that three scaling regimes might be discerned Q. 

At least, the value of the transition scale from in- 
homogeneity to homogeneity needs to be tested 0, . 

We use the more recent and complete three- 
dimensional galaxy catalog, the Sloan Digital Sky 
Survey (SDSS), to repeat older calculations and hope 
to obtain more reliable results. 

We leave out of the account some other issues 
related to the galaxy case. These are source evolution 
and cosmological effects |lfj, |llj, ll2( . They should be 
considered in detail in some further study. 



II. THE CORRELATION INTEGRAL 
METHOD 

As a first approach, we have chosen to use a charac- 
terization of the structures of point sets which is given 
by the correlation integral |Hj, denned as: 



C 2 (r) 
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where O is the Heaviside function, 
mation is over the whole set of N 



The inner sum- 
- 1 galaxies with 



coordinates Xj, j ^ i, and the outer summation is 
over a subset of N' galaxies, taken as centers, with 
coordinates Xj. By taking only the inner TV' galaxies 
as centers we allow for the effect of the finiteness of 
the sample Q. 

This characterization is also valid when the set is 
not fractal. Therefore it seems appropriate to use 
this approach to analyse galaxies considered as point 
sets, provided the spaned scale range shows either a 
fractal behavior or not. 

We may interpret C^r) as Af(r)/N where Af(r) is 
the average number of galaxies within a distance r of a 
typical galaxy in the set. As r goes to zero, C 2 should 
vanish as C 2 oc r° 2 . For computational purposes it is 
more convenient to use the form: 



log(C 2 ) = CONST. + D 2 log(r) 



(2) 



The exponent D 2 is the fractal dimension, neces- 
sarily < 3 for an embedding space of dimension three. 
When D 2 is different from three, the distribution is 
fractal Q. 
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FIG. 1: D2 as a function of r, for small scales. 
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FIG. 2: D2 as a function of r, for larger scales. 



III. THE ANALYSIS 

We use the publicly available SDSS data. After 
appropriate eliminations we are left with 84k galaxies. 
The rms galaxy redshift errors are estimated to be 
about 30km/sec [lfjj, therefore they are small for the 
overall density fluctuations. 

We assume Ho — 65 km/s/Mpc. Since our results 
exhibit large error bars, the value of Hq is not too 
important. Another choice would only shift h^ 1 the 
transition scales. 

Because of the large differences between the tran- 
sitions from one dimension to another, we devide the 
appropriate range into a large number of segments 
and calculate D2 by the Procaccia-Grassberger 
method 01 ■ To obtain a very rough estimate on the 
region of transition, we check three points at a while. 
It gives the values for D2 which are shown in Figures 
1 and 2. 

In the above figures, a step variation of the fractal 
dimension D2 is obvious. This could explain why pre- 
vious studies, limited in scale range, concluded to con- 



stant fractal dimensions with different values. How- 
ever, large fluctuations appear at each transition scale. 
Then D2 decreases towards a fluctuating value and 
rises again with large fluctuations at the next transi- 
tion scale. The fluctuations are such that the fractal 
dimension becomes larger than the limiting dimen- 
sion three at each transition. This seems to point out 
to some artefacts or bias due either to the employed 
method or to the data sample. The interpretations we 
propose are twofold. 

A. First interpretation (tentative) 

If one puts aside the low r results of Figure Q 
because of the too unphysical values attained by D2, 
one can use only the results of Figure |21 and proceed 
to their following processing. 

One devides the range from to 200 Mpc into 256 
segments. To obtain a rough estimate of the region 
of transition, one first checks three points at a while. 
This gives 128 values for D 2 (see Figure EJ. The 
error on the average of the so obtained D% is taken 
as the error (more reliable than the (smaller) error 
obtained from the direct least square fit), because the 
distribution of these segment values is more Gaussian 
than the total least square fit. This method tends 
to give large errors on each segments but enables to 
observe more clearly the trend. Then one calculates 
the dimensions by using all the points to perform a 
least square calculation. 

To check these results, we used two different 
amounts of inner points in the evaluation. Still we 
kept in mind that all the points, including those 
outer points which are at a distance used in our 
calculations, should be within the measured range. 
We also used, as the maximal distance considered, 
three different values. In all the cases the transition 
area remained the same. We then took an inner 
region which is much larger so that some of the points 
of the outer region are already in an area not covered 
by the catalogue. For this case the results differed 
considerably. We do not get a clear transition and 
the dimension was on the average 2.3 ± 0.3. The 
trend is obvious as less and less points are counted 
from the total available points set we expect to get 
more and more distorted results. 

In Figure |3| a transition seems to appear between 
two scaling regimes: 
1.90 ± 0.03 , for 6 Mpc to 80 Mpc. 
3.01 ± 0.04 , for 100 Mpc to 200 Mpc. 

Between 80 Mpc and 100 Mpc the picture is not 
clear, so that the transition might be somewhere 
around 90 Mpc. In Figure^] which is obtained by the 
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FIG. 3: The tentative D 2 = 2 to D 2 = 3 transition. Tak- 
ing the parameters of the least square fit for the lower 
part and upper part, compared to the experimental re- 
sults (solid line). 
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FIG. 4: The tentative D2 « 1 to D2 — 2 transition with 
the same procedure as for Figure 3. 



same method applied to the low r results of Figure 

another transition seems to appear between a D2 ~ 1 

and D2 ~ 2 regime: 

0.89 ± 0.01 for 2 Mpc to 4.5 Mpc. 

1.93 ± 0.05 for 6 Mpc to 20 Mpc. 

Between 4.5 Mpc and 6 Mpc the picture is also 
not too clear. At the transition, the dimension in- 
creases very dramatically. This might be interpreted 
as an artefact resulting from the small difference in 
r of neighbouring bins. But these low r analyses can 
also be viewed as some confirmation of the prelim- 
inary results obtained by Q who found, for galaxy 
separations up to about 5 Mpc a D% dimension about 
1.2 and for larger scales (out to about 30 Mpc, which 
was the limiting scale of their study) a dimension of 
about 1.8. 



B. Second interpretation (seemingly more 
robust, but to be confirmed) 

Figures 1 and 2 show a similar behavior, with a 
peak at the transition, then an (e _r )-like decrease 
toward a constant value (with fluctuations above the 
allowed value D2 = 3). When more data become 
available, and with an analysis method more adapted 
to the study of multi-fractal distributions, we might 
be able to check whether this peak is actually an 
artefact coming from our data analysis method (some 
step behavior at each transition scale) or if it is due 
to the data sample. 

The D2 = 3.3 peak at 100 Mpc inland the decrease 
farther could be of the same nature: in this case, the 
"true" value of D2 beyond 100 Mpc would not yet have 
been reached at the limit of the study (200 Mpc) (it 
might be < 2.6), which could mean that the transition 
to uniformity {D2 — 3) has not yet been reached at 
these scales. 



IV. CONCLUSION AND DISCUSSION 

We use the publicly available data from the SDSS, 
to complete an analysis of the fractal dimension of 
the galaxy distribution, with the correlation integral 
method. 

We check scales up to 130 h^ 1 Mpc. We obtain an 
obvious step variation of the fractal dimension D%. 
This could explain why previous studies, limited in 
scale range, concluded to constant fractal dimensions, 
with different values. 

Two possible interpretations of these results are 
proposed: 

• A rough mean square fit method gives i) up to 
4.5 to 6 Mpc, a dimension of the order one, ii) 
then a transition to a dimension of the order 
two, iii) and between 80 and 100 Mpc, another 
transition to a dimension around three. 

• However, the variation of D2 with scales show a 
peak at each transition (the fractal dimensions 
at the peaks become larger than 3), then an 
(e~')-like decrease toward a constant value 
(with fluctuations). This could be due to the 
step behavior at the transition scale. The 
transition to uniformity (D2 — 3) would thus 
not yet have been reached at the largest studied 
scales. Moreover, a transition to homogeneity 
at 130 h^ 1 Mpc would be inconsistent with 
the sizes of the largest structures seen in the 
universe 0, 0|. Last, it would be interest- 
ing to check if a variable with scale fractal 



4 



22nd Texas Symposium on Relativistic Astrophysics, Stanford University, December 13-17 2004 



dimension of the galaxy distribution might 
be related to theoretical predictions proposed 
as a consequence of a principle of relativity 
of scales 18]: a transition to homogeneity 
predicted around 750 Mpc and a multi- 
fractal distribution with a dimension varying 
by steps whose width might be related to the 
organization hierarchy observed for the galaxies. 

However, these large fluctuations, which appear 
at each transition scale, seem to point out to some 
artefacts due either to the employed method or to 
the data sample. A way to discriminate between 
the two possible reasons of the appearance of these 
odd fluctuations would be to test the validity of 
the application of the Grassberger and Procaccia's 
correlation integral method to a mutifractal distri- 
bution by applying it to a set of mock catalogues 



of galaxies, artificially constructed, with a known 
fractal dimension varying by known steps. If the 
runs of the here employed code reproduce the known 
features of the distributions, it would suggest that 
the artefacts might be due to the data sample 
analysed here. We would therefore need a better 
sample to go further on and, e. g., apply our method 
to the next SDSS catalogue, presumably publicly 
available in a very near future. If, on the contrary, 
those runs are not able to reproduce the known 
features of the distributions, it would suggest that 
the artefacts might be due to the employed method, 
and we would have to test other ones and find a code 
better fitted to the study of multifractal distributions. 
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